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We have performed time-domain interferometry experi- 
ments with matter waves trapped in an harmonic potential 
above and below the Bose-Einstein phase transition. We in- 
terrogate the atoms according to the method of separated os- 
cillating fields, with a sequence of two radio-frequency pulses, 
separated by a time delay T. We observe the oscillation of the 
population between two internal Zeeman states, as a function 
of the delay T. 

We find a strong depletion of the interference fringes for 
both the Bose condensates and the thermal clouds above con- 
densation, even at very short times, when the clouds are still 
overlapping. Actually, we explain the observed loss of contrast 
in terms of phase patterns imprinted by the relative motion, 
as a consequence of the entanglement between the internal 
and external states of the trapped atoms. 
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The question about the coherence of Bose-Einstein 
condensates (BECs) , and the characterization of their 
phase properties have drawn considerable attention in 
the recent literature. 

The first evidence of a definite phase for weakly inter- 
acting condensates dates back to the experiment of the 
MIT group 1^, where high contrast matter- wave inter- 
ference fringes were observed in the density distribution 
of two freely expanding BEC's. Subsequent experiments 
performed at JILA measured the relative phase of 
two condensate in different internal (hypcrfine) states ex- 
periencing almost the same trapping potential. 

Other experiments have further investigated this sub- 
ject [^m^,^. In particular, the group of W. Phillips has 
recently measured the evolution of the spatial profile of 
the phase of BEC's, by using a Bragg-interferometer [Q. 
They have shown that a trapped condensate is character- 
ized by a uniform phase, and after release from the trap 
develops a non-uniform phase profile. 

In a recent Letter |^ we have demonstrated an exper- 
imental method for a sensitive and precise investigation 
of the interaction between two Bose-Einstein condensates 
pO| . In this paper we present an interferometry experi- 
ment performed in the time-domain on the same system, 
which allows to gain a deeper insight into Ramsey in- 
terferometry with ultra-cold atoms across the BEC 
phase transition. 

It is well known that Ramsey fringes are readily observ- 
able with a sample at room temperature. The reason is 
that Ramsey signals rely on the persistence of coherence 



between two distinct atomic levels, no matter what the 
state of the particles center-of-mass is. This is true only 
when the internal and external degrees of freedom are un- 
coupled. Whenever the latter condition fails, a depletion 
of the Ramsey fringes visibility can occur and has indeed 
been observed Thus, the entanglement of external 
and internal states of atoms trapped in a magnetic po- 
tential is the basis for the use of Ramsey method to char- 
acterize the phase properties of a condensate and a ther- 
mal cloud near BEC. In the JILA experiment jsj, such an 
entanglement caused the observed Ramsey fringes in the 
population of the two hyperfine levels to undergo a loss of 
contrast, as a consequence of the reduced spatial overlap 
between the two BEC's, due to their mutual repulsion. 

In our experiment we assist to a similar depletion but 
on a much faster time-scale, of the order of tens of fis, 
that cannot be explained with the reduction of spatial 
overlap. We prepare, and subsequently probe the sys- 
tem with a sequence of two identical radio-frequency 
(rf) pulses, separated by a delay time T. The relative 
phase accumulated by the condensates produces Ramsey 
fringes in the population of each level, as a function of the 
delay T. We observe a strong reduction of the oscillation 
amplitude even at very short times, when the condensates 
are still almost completely overlapping, due to the rela- 
tive velocity acquired by the condensates. We compare 
this reduction with that of a cloud of thermal atoms, at a 
temperature three times larger than the phase-transition 
temperature Q ~ 30c. 

We prepare a condensate of typically 2 x lO'"^ *^Rb 
atoms in the \F = 2,mf = 2) hyperfine level, confined in 
a 4-coils loffe-Pritchard trap elongated along the z sym- 
metry axis fl^ . The axial and radial frequencies for the 
|2) state are ujz2 = 27r x 13 Hz and uj±2 = 2ttx 131 Hz, re- 
spectively, with a magnetic field minimum of 2.86 Gauss. 
Then, we apply an rf-pulsc to split the initial condensate 
into a coherent superposition of different Zeeman |to/) 
sublevels of the F = 2 state. The atoms transferred in 
a different sublevel move away from the |2) equilibrium 
position with an acceleration that depends on their rrif 
value. Thus, the state of motion becomes entangled with 
the internal atomic state. 

In this experiment we use a 24 cycles rf-pulse at 2 MHz, 
which quickly leaves the |2) state with ~ 41% of the 
initial number of atoms, transferring an equal part of 
atoms 41%) to the |1) state and ^ 15% of the atoms 
the |0) state, - 3% to the | - 1) and ~ 0% to the | - 2) 
states. 

Even though all the five levels should be taken into 
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account in the early stages of the evolution (before the 
atoms in the |0), | — 1) and | — 2) states leave the trap), 
the basic features can be explained by considering only 
the dynamics of the two most populated levels, i.e. |1) 
and 1 2). We will discuss later how the other three levels 
affect the overall behaviour of the system. 

We describe our double condensate system by a spinor 
wavefunction, where the upper and lower components 
refers to states |2) and respectively. Immediately 
after the first rf-pulse, the wavefunction can be written 
as: 

o).^^^fia(l) (1) 

where \/Nou{r) is the equilibrium wavefunction of the A^o 
atoms of the initial |2) condensate, in the Thomas-Fermi 
(TF) regime . As the rf-pulse is much shorter than the 
oscillation periods of the harmonic trap, we could safely 
take the spatial wavefunction u(r) to remain unchanged 
and flop only the internal state. 

In the subsequent free evolution, the relative phase be- 
tween the two components accumulates with a rate pro- 
portional to the difference of chemical potentials fi2 — fJ-i- 
Moreover, u{r; 0) is no longer the equilibrium wavefunc- 
tion for 1 2) nor for |1): the spatial wavefunctions evolve 
as dictated by two coupled GP equations into u{r;t) 
and v{r;t), respectively. 

By applying a second rf-pulse, identical to the first, af- 
ter a time delay T, we suddenly mix the two components 

l'^'^^^)- 2 [iu{r;T)e-'^°^ + v{r;T) J 

with LUo = (/i2 - ^ii)/fi. 

We then separate the two internal states and indepen- 
dently count the number of atoms in each, given by the 
integrated square modulus of the corresponding compo- 
nents, showing an oscillating behaviour at the frequency 

For times T much shorter than the harmonic peri- 
ods, we can neglect all but the most relevant effect of 
the spatial wavefunctions evolution: due the differential 
gravitational "sagging", the |1) condensate acquires a 
downward time-dependent momentum —hq{t). In par- 
ticular, we will also neglect the loss of spatial overlap 
arising from the relative displacement. Then, taking 
w(r;T) = u(r;0), v{r;T) = iu{r]0) exp{-iq{T)y) and 
given the normalization J \u{r; 0)\^dy = 1, we have 

N2{T) ^ No^[l ~ A, (T) cos{looT)] 

Ni{T) ^ No ~ N2{T) (3) 

with the slow time-dependent oscillation amplitude 

A,{T)^ I |u(r;0)pcos(g(T)2/)dr. (4) 



The amplitude of the quadrature component vanishes be- 
cause |u(r; 0)p is an even function of the y coordinate. 

The amplitude Ac{T) decays as the relative velocity 
increases and becomes of the order of the vertical 
extension of the original condensate. Eventually, Ram- 
sey fringes are completely washed out. The relative dis- 
placement would give the same result, but only at later 
times. 

In Fig. |l| we plot the experimental data of the oscilla- 
tion amplitudes for the fraction N2/{Ni + N2) versus the 
time delay T. To obtain each data point, we have sam- 
pled one oscillation period around the reported T value. 
We have taken about 10 points per period and each point 
is the average of few acquisitions. To compare the exper- 
imental points with the values predicted by Eq. (Q), we 
rescale the latter by a suitable factor, chosen to match 
the data around T = 0. As for the damping time, we find 
a satisfactory agreement between theory and experiment. 

Actually, for the double condensate we can refine the 
above model to include the mean-field repulsion in the 
wavefunctions evolution: to this end, we numerically in- 
tegrate two coupled Gross-Pitaevskii (GP) equations, ac- 
cording to the model described in Ref. ||T^. On general 
ground, in the first stages of the evolution the phase of 
each component can be written as a quadratic form in 
the spatial coordinates |@,|l5|, describing the mean- field 
expansion/contraction of the condensates in the Thomas- 
Fermi regime |l5| ] and the motion along the vertical di- 
rection y. The results of the numerical simulations, in 
the time range considered here, show the following: (i) 
the phase of the |2) condensate remains almost uniform, 
and only for later times {t ~ 0.5 ms) develops a negative 
curvature, due to the fact that the condensate is con- 
tracting after the initial transfer of atoms to the other 
levels; (ii) the phase of the |1) condensate is dominated 
by a linear term 

771 

Mr,t)-j{vo{t)+6v{t))y (5) 

where vq (t) is the velocity acquired during the fall in the 
trapping potential 

and 6v{t), which is negligible for t < 0.5 ms, is due to the 
mutual repulsion between the two condensates. Thus, the 
GP simulations confirm that, at short times, the simple 
model adopted above correctly describes the basic fea- 
tures of the system. 

The phase term (^ is responsible for washing out com- 
pletely the Ramsey fringes even at very short times, when 
the condensate are still almost completely overlapping (as 
shown in Fig. ^ the overlap between the two wavepackets 
is substantial even at T as long as 0.5 ms). Physically, 
across the spatial extension of the wavepacket there are 
regions of alternated positive and negative interference, 
with a vanishing net transfer of atoms. 
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For a thermal cloud we can consider the system to be in 
a given quantum state and then take an ensemble average 
over all the accessible states. This way, we need only to 
replace the time-dependent amplitude with 

Ath(r) - l^E / ■ l*{«}(r-)pcos(g(T)y)dr (7) 

where /{-„} = [exp((e{-„j — /i)/fc9) — 1]^^ is the Bose 
mean occupation number of the harmonic trap eigen- 
state = ipndx)'^n2iy)ipn3{z) with energy ^ 

is the chemical potential, given by the normalization 
S{n} /{"} ~ -^0, and Q the temperature. By carry- 
ing out the integration over x, z and the sum over the 
corresponding quantum numbers ni, 713, we find 



Ah(r) 



X / \i)nM? cos{q{T)y)dy, 



))x 



(8) 



the g2{x) = X^z^i ^V^^ function being the result of re- 
placing the discrete sum over 711,713 with a double inte- 
gral. In principle, one should allow for the rf-pulses to 
act differently on the different harmonic oscillator eigen- 
states as the detuning varies. This would introduce {71}- 
dependent weights in Eq. (^). However, we have verified 
that for the relevant levels at 8 ~ 0.4 /iK, these weights 
are equal within a few percent. 

We note here that Eq. d7h can be rewritten 



^th(T) = n(r)cos(g(T)y)dr 



(9) 



where n{r) is the spatial density. In this form, it is 
evident that the thermal cloud behaves as a coherent 
wavepacket, the reason being that the displacement oc- 
curring between the two rf-pulses is much less than the 
thermal coherence length (Ath = h/y^ZTTmkQ). 

We point out that for the non-condensed sample, where 
typical densities are a factor 30 lower than those of BEC, 
we neglect the atom-atom interactions, which justifies the 
above single-particle analysis. In Fig. |^ we plot the am- 
plitude of the observed ^^2/(^2 + -^1) oscillation for sev- 
eral values of the time delay T between the two rf-pulses 
and we compare with the corresponding predictions given 
by Eq. (0), rescaled to match the experimental values for 
T ~* 0. As for the damping time, we have again a satis- 
factory agreement; thus, we believe that the cause of the 
Ramsey fringes loss of contrast is well understood. 

However, the experimental values are about a factor 
2 less than those given by Eq. (^. As we already no- 
ticed, part of the atoms end up in the ruf — 0,-1,-2 
Zeeman sublevels after the second rf-pulse. Two are the 
main consequences: (z) the peak-to-peak oscillation am- 
plitude of the N2/{Ni + N2) fraction cannot exceed 0.87, 
both for the condensates and the thermal clouds; (ii) 



the oscillation waveforms deviate from a pure sinusoidal 
behaviour. 

To conclude, we point out that our system is suit- 
able to study the possibility of a revival of the Ram- 
sey fringes after the condensates have been spatially well 
separated. As the center-of-mass of the |1) condensate 
undergoes harmonic oscillations around its equilibrium 
position P,[l0|, it comes back to rest at its initial po- 
sition. According to the above description one should 
expect a revival of the oscillations in the relative popu- 
lation when the condensate |1) and |2) come to overlap 
again with almost vanishing velocity. The numerical so- 
lution of the GP equations of the two-level model shows 
indeed that over time scales of tens of ms the Ramsey 
fringes are characterized by collapse and revival |l7| . 

To summarize, we have studied the Ramsey interfer- 
ence of atomic clouds across the BEC phase-transition, 
with a system where the two involved states have equilib- 
rium positions far apart. We have shown that the phase 
pattern imprinted on the moving wavepacket by its ac- 
quired velocity washes out the Ramsey oscillations of the 
fractional populations well before that the spatial overlap 
decreases. By repeating the same experiment on a ther- 
mal cloud at three times the condensation temperature, 
we have observed that the same mechanism is responsible 
for an even faster damping of oscillation contrast. In this 
respect, there is no substantial difference from a thermal 
cloud and a Bose condensate. The analogy with light 
optics is straightforward: indeed, to observe interference 
between two paths we need only the difference between 
the paths length not to exceed the coherence length. 
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We note that, unlike in usual Ramsey experiments where 
the signal oscillates at the frequency difference between the 
atomic frequency and the interrogation frequency, here we 
record oscillations at the atomic frequency. The reason is 
that the phase of the second pulse remains the same, no 
matter what is the time delay T. Thus, varying T we just 
scan the phase of the atomic dipole. 

We defer a more complete discussion on this point to a 
forthcoming publication. 




Time (jxs) 

FIG. 1. Peak-to-peak amplitude of the oscillating observ- 
able N2/{N2-\-Ni), i.e. the |2) condensate population normal- 
ized to the total number of atoms in |2) and |1): experiment 
and calculation (solid line) 
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FIG. 2. Peak-to-peak amplitude of the oscillating observ- 
able N2/{N2 + Ni), for thermal atoms at 6 = 0.4 ~ 3ec: 
experiment and calculation (solid line) 
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FIG. 3. Overlap of the two condensates at the second 
rf-pulse. Density profiles along the vertical axis y of conden- 
sates |2) (dashed-line) and |1) (solid line) before the pulse, 
and of condensate |1) (dotted line) immediately after. The 
curves are obtained by solving the GP equations for the 
two-state model in Ref. [^. Lengths are given in units of 
a±2 = [h/{muj_i_2)]^^'^ = 0.94 /im; T = 0.5 ms. 
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